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A signature of the internal reionisation of the Milky Way? 
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ABSTRACT 

We present a new semi-analytical model of the population of satellite galaxies of the 
Milky Way, aimed at estimating the effect of the spatial structure of reionisation at 
galaxy scale on the properties of the satellites. In this model reionisation can be either: 
(A) externally-driven and uniform, or (B) internally-driven, by the most massive pro- 
genitor of the Milky- Way. In the latter scenario the propagation of the ionisation front 
and photon dilution introduce a delay in the photo-evaporation of the outer satellites' 
gas with respect to the inner satellites. As a consequence, outer satellites experience a 
longer period of star formation than those in the inner halo. We use simple models to 
account for star formation, the propagation of the ionisation front, photo-evaporation 
and observational biases. Both scenarios yield a model satellite population at z=0 
that matches the observed luminosity function and mass-to-light ratios. However, the 
predicted population for scenario (B) is significantly more extended spatially than for 
scenario (A), by about 0.3 dex in distance, resulting in a much better match to the 
observations. The survival of this structural signature imprinted by the local UV field 
during reionisation on the radial distribution of satellites makes it a promising tool for 
studying the reionisation epoch at galaxy scale in the Milky Way and nearby galaxies 
resolved in stars with forthcoming large surveys. However, more work is needed to de- 
termine how the effect reported here can be disentangled from that of cosmic variance 
between different realisations of Milky Way haloes. 
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- cosmology: reionization - diffuse radiation 



INTRODUCTION 

Reionisation took place as the first stars of the universe 
brought the dark ages to an end. It seems now well es- 
tablished that reionisation also had a profound impact on 
galaxy formation. The metagalactic UV radiation field may 
be responsible for evaporating the gas of low-mass galax- 
ies (|Gnedinll2000l : lHoeft et alj|2006h . affecting their star for- 
mation and therefore also the buildup of the galactic halo 
<|Bekki fc Chiball2005l ). This seems to provide an acceptable 
solut i on for the "missin g satellites problem" l|Klvpin et al.l 
1 19991 : [Moore et alii 19991 ). by i nhibiting star forma tion in low 
mass galaxies at early times l|Bullock et al.|[2000l ). 

Several recent semi-analytical models (SAMs) have 
shown that the observed luminosity function (hereinafter 
LF) of the satellites of the Milky Way (hereinafter 



120091: IWalker et all 120101 ). They suggest that the ultra- 
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20081), although the latter is hotly debated ( Aden et al 



Zucker et aimOOfj; iBelokurov fc et al.l 120071 ; llrwin fc et al l 
120071 ; IWalsh et alj 12007) are reionisation fossils, living in 
sub-haloes of about 1O 6_9 M . The success of these mod- 
els invites us to reconsider UFDs, not as a problem, but as 
probes of the reionisation epoch. 

Most past studies implemented a uniform, instan- 
taneous reionisation by simply sup pressing star forma- 
tion in the sub-haloes below z re ion (i Koposov et al.l 120091 : 
iBusha et ail l20ld : iMaccio et aill2010l ). Others tracked the 
evolution of the physical properties of the gas (temperature, 
neutral fraction) subject to the UV background from exter- 
nal and local sources, and SuperNova (SN) feedback, in or- 
der to d etermine self-consistently the satellit e 's reionisation 
redshi ft (|Munoz et al.ll2009l : lFont et al.ll2oTll) . lLunnan et al.1 
l|201ll) studied the effect of a complex, patchy reionisation 
on the satellites properties. In this letter, we investigate the 
effect of varying the structure of the UV field at reionisation 
on the properties of the satellites of the MW. We choose 
2 deliberately simple configurations: a uniform background, 
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and a unique, galacto-centric source. We show that the latter 
produces a significantly more extended satellite population. 

We first present our model and its ingredients (Sec. [2]), 
and then our results in terms of LF, masses, and radial dis- 
tribution of satellites (Sec.|3|. Our conclusions are presented 
in Sec. H 



A SEMI-ANALYTICAL MODEL OF THE 
SATELLITES OF THE MILKY WAY 



Our SAM is partly inspired by iBusha et"al] |20ld ) and 
iKoposov et~aiT (|2009l ) (hereinafter BIO and K09), but de- 
scribes the reionisation epoch in more detail. We allow for 
a non-steady, non-uniform UV background and compute in- 
dividual evaporation times for the gas hosted by the dark 
matter sub-haloes. 



2.1 Ingredients 

The life of a model dark matter sub-halo in the range 
1O 7-9 M0 (typical of a MW satellite at z=0) can be summa- 
rized as follows: after its formation, it grows until it reaches a 
mass sufficient for forming stars. During reionisation, ionisa- 
tion fronts (hereinafter I-front) of galactic or external origin 
propagate, and eventually overtake the sub-halo at a redshift 
zi. The latter starts to lose its gas via photo-evaporation, 
which rapidly suppresses star formation, between 0.4-2 Gyr 
after the Big Bang. Then it undergoes a purely passive evo- 
lution for the next « 6 Gyr, until it is accreted by the MW, 
around z ~ 1. Depending on its impact parameters and 
concentration it will be tidally and ram-pressure stripped. 
Finally, it may get disrupted, losing all of its stars to the 
halo, or conserve a fraction of its pre-accretion mass and 
stellar content, and survive as a satellite until today. The 
purpose of the SAM described here is to model this chain of 
physical processes, with a number of deliberate simplifica- 
tions, in order to produce a population of satellites and its 
observables. The major ingredients are listed below, in the 
order of the modelling steps. 

(i) Sub-halo tracks: a s in BIO we use the Via Lactea II 
(hereinafter VLII) data (|Diemand et al.|[2007l '). The simula- 
tion represents a MW-size halo (about 2.10 12 Mq) that has 
not had any recent major merger, making it a suitable host 
for a MW-like disk galaxy. The dark matter particle mass 
is 4.10 3 Mq, which allows about 20000 sub-haloes to be re- 
solved, down to a few times 1O 5 M0. 

(ii) Threshold mass for star formation M t : a sub-halo is 
allowed to form stars when it reaches T v i r = 8.10 3 K, which 
translates to M t « 1O 7_8 M0. This reflects the typical halo 
mass where atomic hydrogen cooling becomes efficient, as in 
B10. 

(iii) Following B10, once the threshold mass is reached, 
the halo forms stars at a rate SFR = e(f C oidM7) a , with e = 
10" 5 , fcoid = fib/fi m , a = 2, and M 7 = M/1O 7 M , where M 
is the dark matter mass of the sub-halo. 

(iv) Reionisation redshifts: we define the reionisation red- 
shift zi as the time where the I-front reaches the sub-halo. 
The details for both external and internal scenarios are sum- 
marized in Sec. 12.21 At variance with most of the literature, 
star formation does not end at zi, as photo-evaporation is 
not considered instantaneous here. 



(v) Photo-evaporation: we compute evaporatio n time s for 
the gas in the haloes using Eq. (5) of llliev et al.l (|2005l ). as 
a function of the normalized flux Fo, the dark matter mass 
M of the halo, and its collapse redshift z co ii: 



tev = AM7 Fq +D l0S 



o(Fo) 



1 + z c , 

To~ 



Myr, (1) 



where (A,B,C,D,E,F) are taken from llliev et all (|2005l ) for 
the case of a 5.10 4 K black body. Expectedly, t cv increases 
with increasing mass and decreasing UV flux. For simplic- 
ity we assume a constant z co n = 15 for all haloes. Photo- 
evaporation starts at zi, and t ov is converted to a Az ov , so 
that evaporation finishes at a redshift z cv = zi + Az ev - The 
SFR of the sub-halo is set to for 2 < z cv . The values taken 
by Fo in our 2 reionisation scenarios are detailed in Sec. 12.21 
(vi) The accretion redshift z acc is computed for each sub- 
halo using the VLII data. We consider that accretion on the 
MW halo results in instantaneous gas removal (if any left) 
from the satellite through ram-pressure stripping. We set 
SFR = for z < z acc . We can thus write the SFR(M /yr) 
prescription more synthetically as: 



SFR: 



e(fcoidM 7 ) Q if 
otherwise . 



T vir ^ 8.10^ and 

z > max(z cv , z acc ) (2) 



This formalism gives rise to typical periods of star formation 
of several 100 Myr between z = 15 - 6, with SFR=10" 6 - 
lO _2 M0/yr for the observable satellites. 

(vii) Tidal stripping and disruption: A halo is considered 
disrupted if it loses more than 99% of its pre-accretion mass. 
However, for a surviving halo we assume that the stellar 
content re mains untouched. Th is is similar to the approach 
of K09 and lMaccio et all ^OloD . 

(viii) Stellar population models: integrating eq. @ over 
time yields the stellar mass and the individual star forma- 
tion hist ory (SFH) of the sub - haloes . We use the latter along 
with the iBruzual fc Charlotl ([20031 ) models to assign lumi- 
nosities to our satellite population, assuming a constant low 
metallicity Z=0.0004. 



2.2 Reionisation scenarios 

We consider two scenarios, where reionisation of the IGM 
and the subsequent evaporation of the satellites' gas is 
driven either by an external or a galacto-centric source. Note 
that here the 2 scenarios are mutually exclusive. 



Scenario (A): external reionisation 

In this scenario, the UV background responsible for photo- 
evaporating the gas of the satellites is of extragalactic origin. 
We assume it is uniform, as in a majority of previous studies, 
and appears at zi = 15 for all sub-haloes, corr esponding to 
the b eginning of the reionisation of the IGM (|Kogut et al.l 
2003). We tune the value of the background flux to a con- 
stant Fo = 0.03, in the units of llliev et al.l l|2005h . so as to 
match the observed LF of the satellites. This translates into 
J21 = 0.03, which i s comp atible with the fluxes predicted by 
lAubert fc Tevssierl (|2010h for the considered redshift range, 
although in the upper envelope. 
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Scenario (B): internal reionisation by a 
galacto-centric source 

While reionisation app ears patchy in large scale computa- 
tions (Zahn ct al. 20Q3), the morphology of I-fronts at the 
scale of the MW progenitor (5 comoving Mpc) is rather un- 
clear, and de pends strongly on the Lyman continuum escape 
fraction f esc (jlliev et al.ll201l |), which is currently very un- 
certain. Here we assume that star formation within the most 
massive progenitor of the MW provides the photons respon- 
sible for evaporating the satellites' gas and the intervening 
IGM. This model results in a radially decreasing reionisation 
redshift profile i n the IGM around t he emitter, similar to 
those of Fig. 1 of lLunnan et alJ (|201lh . However this model 
can not account for any stochasticity of the reionisation his- 
tory of the IGM. 

(1) Initial I-front propagation: we mo del the I-front 
Ri(z) as a cosmological Stromgren sphere ()Barkana fc Loebl 
l200ll ) with a clumping factor C=l. The SFR of the 
gal acto-centric sou rce is a n extrapolation from Fig. 4 
of ISalvadori et (|201Ch : SFR(z) = SFR 25 10 (25 - z)/10 , 
with SFR25 = SFR(z = 25) = O.O2M /yr. We assume 
that a rate S FR=lMp)/yr produce s 1.6 10 53 ionising 
photons/s as in IShapiro et all (|2004h . We set the start 
of star formation at z — 25, as suggested by theoretical 
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studies llBromm et al. 20091). Nume rical simulations from 
iRazoumov fc Sommer-Lars^l i|2010h show that at z> 10, 
f csc is close to 1. Therefore we use a constant f csc = 1 
throughout the model. 

(2) Photo-evaporation of the sub-haloes: for a satellite at 
a distance Rgc from the Galactic center, photo-evaporation 
starts at zi, when the expanding I-front reaches the satel- 
lite, i.e. Ri(zi) = Rgc(zi)- We convert our SFR to Fo using 
Fp = 16 SFR(z) (10 kpc/RGc) 2 i following the conversions of 
IShapiro et ail l|2004h . 



2.3 Observational sample and biases 

We consider the same observed satellite population as in 
Tab. 2 of K09, consisting of the 11 classical and the 10 re- 
cently discovered SPS S satellites. T h e gala cto-centric dis- 
tances are taken from [Kroupa et all {2010). The observed 
luminosity function is corrected for t he sky coverage of DR 5 
and the detection efficiency given bv lKoposov et all {2008). 
Detected model satellites are required to be within 260 kpc 
of the galactic center (the maximum distance where the tip 
of the red giant branch is detected by SDSS). Moreover they 
need to be within the completeness radius at a given magni- 
tude as given by eq.(8) of K09 at z=0. Finally, their detectio n 
efficiency is computed using eq.(8) of iKoposov et all IMPOST ). 



3 RESULTS 

The LF is degenerated with respect to the parameters 
(e, Fo): a decrease in star formation efficiency e can be com- 
pensated by a longer star formation period, allowed by a 
smaller UV background flux. This is similar to the e — z re 
degeneracy reported in B10. In the present paper, Fo was set 
so as to achieve evaporation of the satellites at high redshift 
on average, i.e. < z cv >~ 6. 
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External 
Internal 
Observed — 
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Figure 1. Observed and model cumulative luminosity functions. 
The shaded areas show the observational error bars corresponding 
to Poisson noise, corrected for biases for the SDSS satellites. Both 
models with reionisation provide a good fit to the observations. 



3.1 Luminosity function 

Model and observed cumulative LFs are compared in Fig. 
[1] Both models with reionisation agree equally well the ob- 
servations. At he faint end, they differ slightly, but are both 
equally likely given the available data. At the bright end 
we underestimate the luminosity of the satellites, because 
of missing physics for the high-mass objects. For instance, 
MW's most massive satellites, such as the large and small 
Magellanic clouds, are still forming stars at z=0, which is for- 
bidden in our model due to the shutdown of star formation 
subsequent to accretion on the MW halo and ram-pressure 
stripping. Our internal and external reionisation scenarios 
are indistinguishable when considering the existing data on 
the LF alone. To illustrate the overall effect of reionisation, 
we also show the results of a model identical to model (A) 
but without reionisation. It overestimates the abundance of 
satellites by a factor 10 for all but the 2 brightest objects 
and is th erefore not favo red in our framework. We refer the 
reader to lLi et all (|2010l ) for an in-depth exploration of mod- 
els without reionisation. 



3.2 Masses 

The central masses (M300, within the 300 central pc) of the 
model satellites are compared with observations in the top 
panel of Fig. [2] They are taken directly from the VLII data 
at z=0. Both reionisation scenarios r eproduce the range o f 
values of M300 and My derived by IStrigari et alJ (2008). 
The distribution of M300 is almost flat, but there is a 
slight M300 — My correlation which is not clearly seen in 
the data. This is typical of this cl ass of models and i s 
seen in all published works so far ilMufioz et al.l {2009); 
iMaccio et ail (|201Cj ): iFont et all (|201lh :K09:B10'). as noted 
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Figure 2. Masses of the faint satellites. Top: central masses 
(M300), models and as derived by Strigari et al. (2008). Bottom: 
total masses of simulated satellites. 

bv lKroupa et all l|2010l ). Eventhoug h the physical mea ning 
of this relation is currently unclear (I Walker et al.ll201(ih . we 
include it as a comparison with previously published mod- 
els. The bottom panel of Fig. [2] shows the total dark matter 
masses of the satellites: despite the apparent common mass 
scale of M300 « 1O 7 M0, the satellites do indeed span a wide 
range in total mass, with a large scatter. As for the LF, the 
masses of the satellites do not allow discrimination between 
external and internal reionisation. 



3.3 Radial distribution of satellites 

In order to quantify the impact of the limited sky coverage 
of the SDSS (about 1/5 in area), we produced mock sur- 
veys using our model satellite population: we include the 11 
brightest to simulate the "classical" satellites and randomly 
pick 10 of the fainter ones within the detectable subsample 
(i.e. the ones that make up the LFs of Fig. [TJ). We show 
the mean (solid lines) and dispersion (shaded areas) of the 
cumulative radial distributions at z=0 obtained for 1000 of 
these mock surveys in Fig. [3] There is a strong difference in 
the predicted profiles of our two reionisation scenarios: the 
distribution for the internal reionisation model is shifted out- 
wards by about 0.3 dex with respect to the external reioni- 
sation model. This shift is 2 times larger than the dispersion 
of the predicted profiles, and is therefore significant. This is 
the signature of internal reionisation: in the galacto-centric 
reionisation scenario the satellites of the inner halo see a 
more intense UV flux than their outer halo cousins, and 
therefore evaporate faster. Thus the outer halo satellites ex- 
perience a longer star formation activity period and end up 
brighter. In this model, all the detected satellites beyond 
R>150 kpc see their photo-evaporation achieved later than 




log(RGc[kpc]) 

Figure 3. Cumulative radial distribution profiles, observed 
(black), and predicted by the internal and external reionisation 
models at z=0. The shaded areas show the dispersion of profiles 
obtained by 1000 realisations of a mock DR5 survey. For better 
readability, the blue area is plotted on top of the red one, thereby 
partially masking it at the top and bottom of the figure. 

z=6, the epoch o f end of reionisation of the IGM derived by 
iFan et al l l|2006l ). 



4 CONCLUSIONS AND DISCUSSION 

We developed a SAM of the satellite population of the MW 
based on the VLII simulation, simple recipes for star for- 
mation and radiative effects, and taking into account the 
observational biases. Our treatment of the propagation of 
the I-front within the early MW progenitor and photo- 
evaporation allows us to explore the effect of the struc- 
ture of the UV field during reionisation on the properties of 
the satellite population. In the galacto-centric reionisation 
model, the satellites of the inner halo are exposed earlier to 
a nearer UV source, and therefore evaporate faster than the 
outer halo satellites. As a consequence star formation in the 
inner halo is suppressed sooner, whereas in the outer halo, 
satellites continue forming stars until up to several 100 Myr 
later. The resulting radial distribution of satellites is shifted 
outwards by 0.3 dex in radius with respect to the external 
reionisation model, and provides a very good match to the 
observations. 

The sensitivity of the radial distribution of satellites to 
the structure of the UV field during reionisation makes it 
an interesting probe for the study of reionisation at galaxy 
scale in the MW with Pan-STARRS and upcoming LSST. 
Provided that the VLII halo is representative of the MW, 
our result suggests that star formation at the center of our 
Galaxy is responsible for the photo-evaporation of its own 
satellites and the intervening gas. However this interpreta- 
tion is complicated by cosmic variance in the assembly his- 
tory of the MW, which can also shift the distribution, as 
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shows Fig. 14 of K09. In order to assess the significance of 
our result for the reionisation history of the MW, we will 
need to couple our model to a larger set of high resolution 
simulations of the formation of the MW similar to the V LII, 
such as the AQUARIUS dat a set dSpringel et ail 120081 ). in 
the spirit of iFont et all (|201ll ): iLunnan et al.1 lj201l£ 

Moreover, a number of simplifications were made in this 
work, which need to be improved on. Firstly, the Iliev et al.l 
(2005) photo-evaporation simulations were designed to be 
mostly reliable around 1O 7 M0. Thus, we need to check via 
further numerical simulations that the evaporation times we 
extrapolated up to the > 1O 8 M0 regime remain realistic. 
We will also investigate the effect of endogenous reionisation 
(i.e. induced by stars formed within the dwarf galaxy itself), 
which has been neglected in this work. Secondly, we need 
to generalize the present framework to account for multiple 
sources in the internal scenario: some fraction of the less 
massive progenitor haloes may also contribute to the UV 
field. Even though they are clustered around the most mas- 
sive progenitor, they could lead to a more patc hy UV field 
aroun d the proto-MW, as seen at larger scale in lZahn et al.l 
The transient nature of the radiation field resulting 
from the short lifetime of massive stars is another untested 
aspect of the model. Accurate evaluation and modelling of 
these effects is crucial for the interpretation of Pan-STARRS 
and LSST surveys of the satellite galaxies of the MW. Since 
the latter are major pro genitors of MW ste llar halo and 
streams, as illustrated in iGomez et aD (|2010l ). accurate in- 
terpretation of Gaia data will render this detailed under- 
standing of the reionisation of dwarf galaxies mandatory. 
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